Thermoelectric Modeling of the Non-Ohmic Differential Conductance in a Tunnel 

Junction containing a Pinhole 
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To test the quality of a tunnel junction, one sometimes fits the bias-dependent differential conduc- 
tance to a theoretical model, such as Simmons's formula. Recent experimental work by Akerman 
and collaborators, however, has demonstrated that a good fit does not necessarily imply a good 
junction. Modeling the electrical and thermal properties of a tunnel junction containing a pinhole, 
we extract an effective barrier height and effective barrier width even when as much as 88% of the 
current flows through the pinhole short rather than tunneling. A good fit of differential conductance 
to a tunneling form therefore cannot rule out pinhole defects in normal-metal or magnetic tunnel 
junctions. 

PACS numbers: 85.75.Dd, 85.30.Mn, 73.40.Gk 



I. INTRODUCTION 



With recent rapid progress toward incorporating tun- 
nel junctions into practical devices, such as magnetic sen- 
sors and MRAMji*2i2ii it has become increasingly impor- 
tant to assure junction quality. Pinhole shorts through 
the insulating layer of a nominal tunnel junction con- 
stitute a possible failure mode. As early as the 1960s 
and 1970s, a set of criteria emerged, attributed to Row- 
ell (who also credits Giaever)^ for determining whether 
a particular device was a good tunnel junction. Only 
a few apply when neither metal junction superconducts. 
One is the exponential dependence of the resistance on 
the thickness of the insulating layer over a series of de- 
vices with increasing thickness; however, a simple model 
shows that classical pinholes can mimic this signature of 
tunnelingpi 

Another, more commonly used, criterion is a 
non-Ohmic differential conductance^*^ Simmons^ and 
Brinkman, Dynes, and RowelUS calculated the tunnel- 
ing current I through idealized junctions as functions 
of bias V, the height of the energy barrier, and the 
thickness of the insulator; when differential conductance 
dl/dV shows positive curvature, can be fit to one of 
these forms, and yields "reasonable" values for the bar- 
rier height and thickness, comparable to independent 
measurements, a good junction is sometimes presumed. 
Akerman and collaborators have demonstrated experi- 
mentally that this criterion alone cannot establish the 
quality of a junctioniA*A£ They made a set of junctions 
with one normal and one superconducting lead. Above 
the superconducting-transition temperature, all fit the 
tunneling models well. Below the superconducting tran- 
sition, in some of the junctions the current was sup- 
pressed inside the gap, as would be expected for tun- 
neling from a normal metal through an intact insulator 
into a superconductor. Other junctions, however, instead 
showed the enhancement of current inside the gap typ- 
ical of Andreev reflection, indicating essentially direct 
superconductor-metal contact 



We now demonstrate, in a simple model, how a tun- 
nel junction with a pinhole short might reproduce the 
non-superconducting result, finding that as much as 88% 
of the current might flow through the short and still 
leave enough tunneling to give a good fit to the Simmons 
model. 

A sample known to harbor a pinhole may show a dif- 
ferential conductance with either positive or negative 
curvature^ the latter a heating effect: a pinhole will 
dissipate more heat at greater bias, becoming hotter, 
and, assuming it acts as a metal, conduct less. By con- 
trast, the conductance of a tunneling channel increases 
with bias, as carriers see an effectively narrower trape- 
zoidal barrier The pinhole and the tunneling chan- 
nels, therefore, have opposite effects on the curvature, 
the one tending to make the conductance curve down, 
the other making it curve upward. For parameters based 
on an experimental geometry, we find that the large pos- 
itive curvature of a relatively small tunnel current can 
overwhelm the weakly negative curvature of the conduc- 
tance through the pinhole short, still yielding an excellent 
fit to the Simmons form with apparent barrier width and 
height varying from the true width and height by factors 
of 2-3. 

Although two of the Rowell criteria are found un- 
reliable, the temperature dependence of conductance 
still appears to distinguish good junctions from shorted 
oneg| ii 4 i2 i i8 4 i9 Recently we have proposed an additional 
test, using only electrical measurements at a single tem- 
perature, that may serve both to diagnose the presence 
of a pinhole and to determine its locationiSi 



II. MODEL AND COMPUTATION 

We build a three-dimensional network, each node of 
which is defined by material, position, thickness, and 
length. The voltage and temperature at each node are 
computed iteratively, in turn determining the local elec- 
trical resistivity and thermal conductivity. We extract 
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zero (Kirchhoff 's law), so that for a node A, 



(i) 



where Vj is the voltage at node j (= A, i) and the sum 
runs over the nearest neighbors of A. The resistance B,Ai 
between nodes A and i depends on the average of the 
temperatures Ta and Ti of the two nodes and, in the case 
of the tunneling layer, on the voltage difference Va — Vi 
through Simmons's formulaAS Equation is set as a 
large sparse matrix and solved with the standard SLATEC 
package. Similarly, the total heat flux through a node is 
equal to the heat generated Ohmically: 
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FIG. 1: Normalized differential conductance against bias. 
The dotted line plots differential conductance absent any pin- 
hole; the solid lines trace conductance as the pinhole side 
changes (top to bottom) from 0.1 nm to 1.5 nm in increments 
of 0.1 nm. The separatrix between negative and positive 
curvature corresponds to a pinhole side between 1.0 nm and 
1.1 nm. The bottom-layer temperature is 77K. 



our dimensions from the experimental geometry in ref- 
erence lilt the side of the tunnel junction is 50/mi; the 
bottom electrode is Nb of thickness 80 nm. On top of 
the Nb layer are an Al layer 8.5 nm thick and a 2.0 nm 
A10 x barrier; then Fe (50 nm) is deposited on top. We 
set the barrier height at 0.5 eV. We model the pinhole 
as a metallic (Al) right square prism at the center of the 
tunnel junction. We alternate calculations of the volt- 
age and the temperature at each node until reaching a 
steady state before measuring the effective resistance of 
the junction. 

The electrical boundary conditions fix one edge of the 
bottom metallic layer at ground and one edge of the top 
metallic layer (rotated 90° relative to the grounded edge) 
at some positive voltage (the "input voltage"). The mea- 
sured bias is the difference between the average voltages 
on the two edges opposite the fixed edges. (This sim- 
ulates a four-terminal measurement of the sort we have 
previously modeled in Reference 1201 The measured bias 
is generally somewhat smaller than the input voltage.) 
We hold the bottom metallic layer in equilibrium with a 
fixed heat bath; all heat generated in the junction must 
leave through this surface. All other surfaces are as- 
sumed perfect thermal insulators. (Radiation losses are 
ignored.) 

There are several approaches to non-linear electro- 
thermal modeling. 2 -' 22 Here, we alternate relaxation 
methods for the voltage and the temperature at each 
node. The total electrical current through each node is 



where Ki is the temperature-dependent thermal conduc- 
tivity corresponding to the material of node i, Xai is the 
distance between nodes A and i, and Sai is the cross- 
sectional area spanned by the link between the nodes. 
First, ||TJ is solved; it must be iterated until voltages 
(checked at the pinhole and away from the pinhole) con- 
verge, since Rai through the insulating layer depends on 
the voltage drop. Then the temperature at each node is 
determined through l|2|l. and the procedure is repeated 
until voltages and temperatures both converge. 



III. RESULTS 

We consider runs with the bottom layer held at 77K 
and with the side of the pinhole ranging from 0.1 nm 
to 1.5 nm. For an input voltage of 0.3V, we find the 
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FIG. 2: The effective thickness of the barrier against the side 
of the pinhole. The true barrier thickness is 2.0 nm. The 
bottom-layer temperature is 77K. 
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FIG. 3: The effective height of the barrier against the side of 
the pinhole. The true barrier height is 0.5 eV. The tempera- 
ture is 77K. 

maximum steady-state temperature at the center of the 
pinhole for the largest pinhole to be 86K. This rise of 9K 
is comparable to that estimated in Reference The 
temperature rise rapidly falls off away from the pinhole. 

As illustrated in Figure ^ the differential conductance 
may show positive or negative curvature, depending on 
pinhole size; we graph normalized curves for the pinhole 
side ranging from 0.1 nm to 1.5 nm, with the separatrix 
between 1.0 nm and 1.1 nm. When its curvature is posi- 
tive, we can fit dl/dV to Simmons's formula to extract 
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FIG. 4: Even in the presence of a large pinhole (0.9 nm on 
a side), the curvature of differential conductance can be fit 
to the form of Simmons. The solid line plots the simulation, 
while the dotted line plots Simmons's formula for the best-fit 
barrier height and thickness. The bottom-layer temperature 
is 77K. 



FIG. 5: The percentage of current that flows through the 
pinhole against the side of the pinhole with the bottom layer 
held at 77K (main graph) and 300K (inset). The input volage 
is 0.3V. 



an effective barrier thickness and effective height. These 
effective barrier parameters are shown in Figures [3 and 
13 Without a pinhole, we recover essentially the input 
parameters. As the pinhole grows in size, the fits remain 
quite good (Figure 0}, but the effective barrier parame- 
ters vary: the apparent barrier thickness decreases, while 
the apparent barrier height increases. The decreasing 
apparent thickness has the greatest effect on the over- 
all conductance (which increases), while the increasing 
apparent height tends mostly to flatten the curve of dif- 
ferential conductance. Such trends are consistent with a 
simple WKB treatment of tunneling. 17 We obtain similar 
results with the bottom layer held at 300K; the separa- 
trix between positive and negative curvature occurs for 
a pinhole side between 4.0 nm and 4.1 nm. 

From the above simulations, we obtain the percentage 
of the total current flowing through the pinhole (Fig- 
ure EJ). This percentage will increase with the side of the 
pinhole; the highest value corresponding to a positive 
curvature of dl/dV is 88% at 77K. (The highest value is 
81% at 300K.) 

Our simulations support the contention that a good fit 
of differential conductance to the Simmons form fails, by 
itself, to verify the quality of a tunnel junction. 



Acknowledgments 

We thank Johan Akerman for useful discussions. DAR 
is a Cottrell Scholar of Research Corporation, which has 
funded part of this research. Numerical work was carried 
out at the Research- Oriented Computing Center of the 
University of South Florida. 



4 



* Corresponding author 

1 S. S. P. Parkin, K. P. Roche, M. G. Samant, P. M. Rice, 
R. B. Beyers, R. E. Scheuerlein, E. J. O'Sullivan, S. L. 
Brown, J. Bucchigano, D. W. Abraham, et al., J. Appl. 
Phys. 85, 5828 (1999). 

2 J. M. Daughton, A. V. Pohm, R. T. Fayfield, and C. H. 
Smith, J. Phys. D : Appl. Phys. 32, R169 (1999). 

3 W. Reohr, H. Hbnigschmid, R. Robertazzi, D. Gogl, 

F. Pesavento, S. Lammers, K. Lewis, C. Arndt, Y. Lu, 

H. Viehmann, et al., IEEE Circuits and Devices Magazine 
(Sept.), 17 (2002). 

4 M. Durlam, P. J. Naji, A. Omair, M. DeHerrera, J. Calder, 
J. M. Slaughter, B. N. Engel, N. D. Rizzo, G. Grynkewich, 
B. Butcher, et al., IEEE Journal of Solid-State Circuits 
38, 769 (2003). 

5 J. M. Rowell, in Tunneling Phenomena in Solids, edited by 
E. Burstein and S. Lundqvist (Plenum, New York, 1969), 
p. 273. 

6 D. A. Rabson, B. J. Jonsson-Akerman, A. H. Romero, 
R. Escudero, C. Leighton, S. Kim, and I. K. Schuller, J. 
Appl. Phys. 89, 2786 (2001). 

7 J. S. Moodera and G. Mathon, J. Magn. Magn. Mater. 
200, 248 (1999). 

8 J. Wang, P. P. Freitas, and E. Snoeck, Appl. Phys. Lett. 
79, 4553 (2001). 

9 J. G. Simmons, J. Appl. Phys. 34, 1793 (1963). 

10 W. F. Brinkman, R. C. Dynes, and J. M. Rowell, J. Appl. 
Phys. 41, 1915 (1970). 

11 B. J. Jonsson-Akerman, R. Escudero, C. Leighton, S. Kim, 

I. K. Schuller, and D. A. Rabson, Appl. Phys. Lett. 77, 
1870 (2000). 

12 J. J. Akerman, R. Escudero, C. Leighton, S. Kim, D. A. 
Rabson, R. W. Dave, J. M. Slaughter, and I. K. Schuller, 
J. Magn. Magn. Mater. 240, 86 (2002). 

13 G. E. Blonder, M. Tinkham, and T. Klapwijk, Phys. Rev. 
B 25, 4515 (1982). 

14 G. E. Blonder and M. Tinkham, Phys. Rev. B 27, 112 
(1983). 

15 J. Schmalhorst, H. Briickl, M. Justus, A. Thomas, 

G. Reiss, M. Vieth, G. Gieres, and J. Wecker, J. Appl. 
Phys. 89, 586 (2001). 

16 C. B. Duke, Tunneling in Solids, vol. 10 of Solid-State 
Physics (Academic Press, New York, 1969). 

17 E. Merzbacher, Quantum Mechanics (Wiley, New York, 
1997). 

18 J. J. Akerman, J. M. Slaughter, R. W. Dave, and I. K. 
Schuller, Appl. Phys. Lett. 79, 3104 (2001). 

19 U. Riidiger, R. Calarco, U. May, K. Samm, J. Hauch, 

H. Kittur, M. Sperlich, and G. Giintherodt, J. Appl. Phys. 
89, 7573 (2001). 

20 Z.-S. Zhang and D. A. Rabson (2003), preprint, URL 
cond-mat/0306153 

21 A. L. Palisoc and C. C. Lee, J. Appl. Phys. 64, 6851 (1988). 

22 D. Georg, L. Steffi, and K. Erich, IEEE Transactions on 
very large scale integration (VLSI) systems 5(3), 6515 
(1997). 

23 Simmons gives dl/dV at zero temperature. We have incor- 
porated heuristically an additional linear increase in the 
overall cond ucta nce with temperature following the data 
in Reference I lit a more detailed treatment would take as 
its starting point thermal smearing of the Fermi surface, 



as in J. J. Akerman, I. V. Roschin, J. M. Slaughter, R. W. 
Dave, and I. K. Schuller, Europhys. Lett. 63, 104 (2003). 



